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Anisotropics in the gamma-ray sky from millisecond 
pulsars 



o 

(N 

S3 
33 

oo 

X 
^3 

6 



(N 
> 

o 
in 
in 



Jennifer M. Siegal-Gaskins 1 *, Rebecca Reesman 1 , Vasiliki Pavlidou 2 f, 
Stefano Profumo 3 , and Terry P. Walker 1 

1 Center for Cosmology and Astro- Particle Physics, The Ohio State University, Columbus, OH 43210 USA 

2 Astronomy Department, California Institute of Technology, Pasadena, CA 91125 USA 

3 Department of Physics and Santa Cruz Institute for Particle Physics, University of California, Santa Cruz, CA 95064 USA 



X 



Accepted Received ; in original form 



ABSTRACT 

Pulsars emerge in the Fermi era as a sizable population of gamma-ray sources. Millisec- 
ond pulsars (MSPs) constitute an older subpopulation whose sky distribution extends 
to high Galactic latitudes, and it has been suggested that unresolved members of this 
class may contribute a significant fraction of the measured large-scale isotropic gamma- 
ray background (IGRB). We investigate the possible energy-dependent contribution of 
unresolved MSPs to the anisotropy of the Fermi-measured IGRB. For observationally- 
motivated MSP population models, we show that the preliminary Fermi anisotropy 
measurement places an interesting constraint on the abundance of MSPs in the Galaxy 
and the typical MSP flux, about an order of magnitude stronger than constraints on 
this population derived from the intensity of the IGRB alone. We also examine the 
possibility of a MSP component in the IGRB mimicking a dark matter signal in 
anisotropy-based searches, and conclude that the energy dependence of an anisotropy 
signature would distinguish MSPs from all but very light dark matter candidates. 

Key words: gamma-rays: diffuse background; pulsars: general; methods: statistical 



1 INTRODUCTION 

In the era of precision gamma-ray astronomy, with data of 
unprecedente d quality from the F ermi Large Area Telescope 
' (Fermi-LAT, lAtwood et al J 12009) and ground-based Atmo- 
spheric Cherenkov Telescopes, including H.E.S.S., VERI- 
TAS, and MAGIC, long-standing questions about the high- 
energy universe might soon be successfully addressed. One 
of these is the detailed nature and origin of the diffuse 
gamma-ray emission. The gamma-ray sky is dominated at 
low Galactic latitudes by a bright diffuse Galactic compo- 
nent, stemming dominantly from processes involving cos- 
mic rays such as inelastic hadronic collisions producing neu- 
tral pions, and inverse Compton and bremsstrahlung emis- 
sion from relativistic cosmic -ray electrons and positrons 
(see, e.g., IStrong et af]|200Ch . At high latitudes, the dif- 
fuse gamma-ray background is customarily attributed to 
extragalactic gamma-ray emitters, such as blazars (e.g., 
Istecker fc Salamonlll996r ). Recent Fermi-LAT results, how- 
ever, indicate that resolved blazar s only contribute a small 
fraction of the observed emission ( Abdo et al. 2009b; how- 
ever, see also Abazaiia n et al.l l2010h . in contrast to, e.g., 
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the diffuse X-ray background ([Brandt fc HasingeH 120051 : 
iHickox fc Markevitch|[200^ . l2007D . 

Since the disc overy of period i c ga mma-ray emis- 
sion from pulsars (|Browning et al.l Il97ll ). the possibil- 
ity that this source class contributes non- negligibly to 
the diffuse gamma-ray emission has been considered 



ilBhattacharya fc Srinivasarj Il99ll : iBailes fc Kniffenl Il992l : 
iBhatia et al.l Il997l ). Some of the brightest emitters in the 
Fermi-LAT gamma-ray sky are in fact associated with pul- 
sating objects, often corresp onding to pulsars observed at 
radio and X-ray frequencies (|Abdo et al.ll2009c] ). Compared 
to its predecessor EGRET, Fermi is shedding light not only 
on young, powerful "ordinary" pulsars (with typical rota- 
tion periods of the order of 0.01-1 sec and ages ranging be- 
tween 10 3 and 10 6 yr) but also on a distinct class of periodic 
gamma-ray emitters with much shorter pulsating periods 
(on the order of a few milliseconds), i.e. millisecond pulsars 
(MSPs). The characteristic age r c of MSPs, extrapolated 
from their period and period-derivative, indicates that these 
obj ects are much older than ordinary pulsars, with r c ~ 10 10 
yr l|Abdo et al.ll2009al ). MSPs are thought to be associated 
with binary systems, the spin-up of the pulsar period being 
fueled by accretion of ma ss and angular momentum from the 
neutron star companion l|Phinnev fc Ku lkarni 19941: lLorimerl 
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l200ll ). While the determination of the age of MSPs is a de- 
bated matter given the highly non-trivial nature of thei r 
evolutionary history (see, e.g.. iKiziltan fc Thorsettl [2010). 
the significantly longer lifetime of these objects compared 
to that of ordinary pulsars might offset a birthrate that is 
necessarily lower (given the binary nature of MSPs), and as 
a result the MSP contribution to the Galactic gamma-ray 
luminosity may not be small compared to that of ordinary 
pulsars. 

Despite the dramatic increase in the number of de- 
tected gamma-ray pulsars in the Fermi era, the bulk of the 
pulsar contribution to the gamma-ray sky very likely origi- 
nates fr om a large population of u nreso lved source s. For in - 
stance, iFaucher-Giguere fc Loebl (|201fj| . hereafter 
examined models for the unresolved MSP population and 
found that in some optimistic models the MSP contribution 
to the diffuse background could even be dominant at certain 
energies. In their "viable model" the small set of MSPs de- 
tected by Fermi imply almost 50k unresolved MSPs. The 
gamma-ray emission from ordinary pulsars is very likely 
confined to rather low latitudes (see, e.g.. iHardind Il98ll ; 
iBhattacharva fc Srinivasanl ll99lT ) , reflecting the fact that 
pulsars are born in the Galactic disk, and that ordinary 
pulsars are relatively young objects. On the other hand, the 
product of typical pulsar kick velocities and the character- 
istic age of MSPs implies a length-scale that is much larger 
than the thickness of the Galactic plane, suggesting that 
MSPs should have a broad latitudinal distribution. This is 
reflected in the observed latitudinal distribution of ordinary 
versus m illisecond pulsars d etected by the Fermi-LAT (see 
Fig. 1 in lAbdo et alj|2010ah . In this respect, MSPs can con- 
tribute to the diffuse gamma-ray emission at high latitudes 
where the Galactic diffuse component is generally thought to 
be comparable or sub-dominant with respect to an isotropic 
extragalactic background. 

Interestingly, however, measurements of the spectrum 
of the large-scale isotropic diffuse gamma-ray background 
(hereafter IGRB) by Fermi find that it is consistent with 
a power law at en ergies between 250 MeV and 50 GeV 
|Abdo et al.l l2010bh . while MSP spectra exhi bit a strong 
cut-off feature at typical energies of a few GeV (|Abdo et al.l 
l2009af ). This implies that either the MSP contribution is sub- 
dominant with respect to the primary IGRB component at 
these energies, or that a complicated combination of sev- 
eral components with peculiar spectral features - e.g., a 
star-f orming galaxy com ponent with a feature at ~300 MeV 
(e.g., iFields et al.ll2010l ). a MSP component with a feature 
at a few GeV, and a hard blazar component dominating at 
higher energies - combine in such a way that they appear 
as an overall almost featureless power-law - a contrived sce- 
nario, but one that cannot be excluded in principle. In either 
case, it appears that it will be difficult to detect spectrally 
the presence of a MSP component in the IGRB, although 
it remains possible to put conservative constraints on the 
unresolved MSP ga mma-ray emission based on IGRB mea- 
surements (see, e.g.. IF-GL10u . 

A powerful tool to investigate the nature of diffuse 
emission is to explore the intensity variation of the emis- 
sion in the sky, e.g., via the calculation of an angular 
power spectrum of anisotropies. Recent theoretical work 
has generated predictions for the angular power spec- 
trum of the gamma-ray emission originating from sev- 



eral known and proposed source classes. These include 
confir med extragalactic g amma-ray po p ulatio ns such as 
AGN (lAndo et al] 120071; [Miniati et all 1200 it ) and star- 
forming galaxies (lAndo fc Pavlidoul 12009] ) . as well as dark 
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In addition, since the distribution of dark matter 
subhalos in our Galaxy is quite radially extended, gamma- 
ray emission from annihilation and decay in Galactic sub- 
structure appears remarkably isotropic on large angular 
scales, although the clustering of dark matter in subhalos 
leads to small-scale anisotropies. Consequently, these struc- 
tures may provi de a substantial contr i bution to anisotropies 
in the IGRB (ISiegal-GaskinsI |200S| ; IFornasa et all 120091 ; 
llbarra et aL1l2009l ; lAndoll2009l ). 

The combined use of spectral and anisotropy in- 
formation in the IGRB (the anisotropy energy spec- 
trum) could conceivably help reveal the presence of 
even a subdominant componen t in the diffuse emission 
l|Siegal-Gaskins fc Pavlidoull2009r ). In particular, it has been 
shown that the anisotropy energy spectrum could be a sensi- 
tive pr obe of the presence of a dark matter comp onent in the 
IGRB (|Henslev et al.ll2010l ; ICuoco et al.ll2O10h . This tech- 
nique is also promising for detecting a subdominant MSP 
contribution to the IGRB, since the emission from unre- 
solved MSPs is expected to feature much stronger anisotropy 
than the extragalactic component, due to the fact that MSPs 
are relatively few and nearby, compared to cosmological pop- 
ulations that may constitute the dominant contributors to 
the IGRB intensity. 

Additional motivation to study the gamma-ray 
anisotropy properties of MSPs is provided by the poten- 
tial interference of MSPs with anisotropy-bas ed searches for 
dark matter. Fermi data (|Abdo et al.ll2009al ) indicate that 
the typical gamma-ray MSP spectrum is, in fact, uncom- 
fortably similar in its overall features to what is expected 
for the annihilation or decay of certain particle dark matter 
candidates, especially if the dark matter is light (tubm < 
few tens of GeV). Furthermore, although the amplitude of 
anisotropies from dark matter annihilation is uncertain, in 
some scenarios it is expected to be quite large, and thus 
it is conceivable that a MSP-induced modulation in the 
anisotropy energy spectrum of the IGRB could be confused 
with a similar modulation induced by dark matter. 

In this paper, we explore the potential of an angular 
power spectrum measurement of the IGRB to probe the 
properties of the Galactic MSP population. We demonstrate 
the power of this approach for an example class of MSP 
population models by deriving constraints on those mod- 
els from the Fermi prelim inary aniso tropy measuremen t 
l|Siegal-Gaskins et al. 201p| ; see also I Vargas et al. 201(ih . 
The model prescriptions we adopt to describe the intensity 
and sky distribution of unresolved MSPs are summarized in 
Sj2] and our procedure for generating simulated maps of the 
MSP gamma-ray emission is outlined in In §3] we calcu- 
late the intensity spectrum and energy-dependent angular 
power spectrum of the collective unresolved MSP emission 
for this class of models and discuss those properties in the 
context of other relevant source classes, including dark mat- 
ter. We compare the predicted anisotropy from MSPs to the 
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preliminary Fermi measurement of the angular power spec- 
trum of the IGRB and obtain constraints on the properties 
of the MSP population in iJU We discuss our findings and 
conclude in ij6] 



2 MODELING THE MSP POPULATION 

The properties of the MSP population that affect the mea- 
sured anisotropy are the sky distribution of MSPs and their 
flux distribution. The former is determined by the spatial 
distribution of MSPs in the Galaxy, while the latter is de- 
termined, for a fixed spatial distribution, by the distribu- 
tion of MSP luminosities. In this study we adopt models for 
the gam ma-ray MS P population based on the semi-empirical 
models of lF-GLld . We emphasize, however, that this work is 
a technique demonstration, and therefore its goal is to show 
that MSPs could produce an observable anisotropy signal 
in Fermi-LAT data, and that an anisotropy analysis could 
be used to constrain the collective properties of the Galac- 
tic MSP population; not to perform a detailed study of the 
consistency of a specific model with the data, nor to iden- 
tify which of several models is preferred by the data. With 
that purpose in mind, we fix the values of the parameters 
controlling the spatial and luminosity dist ributions of MSPs 
to those of "viable" model MSP2.base of lF-GLld . and dis- 
cuss the expected impact of variations in these parameters 
on our results in i|6] 

We take the fiducial number of MSPs in the Galaxy 
Nmsp = 49k, as in model MSP2J>ase. Since the observables 
considered in our study (high-latitude intensity and angular 
power) scale straightforwardly with A^msp and the typical 
flux of an individual high-latitude MSP, Pi , we also consider 
the dependenc e of our results on these parameters in £j5] 

Following lF-GL10l . we describe the MSP spatial distri- 
bution with a Gaussian function of radius for the surface 
density projected on the Galactic plane, 

p(r) oc exp(-r 2 /2cr 2 ) < r < 100 kpc, (1) 

where r is the projected distance from the Galactic Centre in 
the Galactic plane, p(r) is the surface density of MSPs, and 
ay, taken to be 5 kpc, characterizes the radial extent of the 
distribution. The latitude distribution of MSPs is assumed 
to follow a simple exponential form, 

N(z) oc exp(-|z|/(|z|}) < z < oc, (2) 

with the scale height (|z|) = 1 kpc. 

Early work on gamma-ray pulsars (e.g., lAronsI Il996l ) 

identified the simple empirical relation L 7 oc between 
the pulsar's gamma-ray luminosity L 7 and the rate it loses 
rotational kinetic energy E — 4ir 2 EP/ P s , where P and P 
are the period and time derivative of the period, respectively, 
and /* is the m oment of inertia of the star. However, recent 
work (see, e.g., lF-GLld ) has found that the luminosities of 
gamma- ray MSPs appear to obey the relation L 7 oc E. As in 
lF-GLld we define the MSP gamma-ray energy luminosity 
(energy per unit time) 

L 7 = min{CP 1/2 p- 3/2 , f?**E}, (3) 

where the proportionality constant C = 10 40,9 erg s 1 ^ 2 and 
jmax _ Q Q5 j g -j- ne assumet [ maximum fraction of rotational 



power loss converted into gamma rays. The integrated pho- 
ton luminosity (photons per unit time) L 7 h above 100 MeV 
is obtained by assuming an energy spectrum with approx- 
imately equal power per decade of energy up to a cutoff 
energy of -B ma x — 3 GeV. For the model adopted in this 
work, it is notable that Eq. [3] results in the vast majority of 
MSPs being assigned luminosities according to the L 7 oc E 
relation. 

For the MSP population, a power-law distribution for 
the rotation period P is assumed, 

N(P) oc P~ 2 1.5 ms < P < 60000 ms, (4) 

and the magnetic field strength B is taken to follow a log 
normal distribution, 

JV(logB) oc exp(-(logS- (log B)f/2a? ofl B ), (5) 

with (log\B) = 8 and a\ OS B ~ 0.2 with B in Gauss. The 
spin-down rate P is determined via the relation B — 3.2 x 
10 19 (PP) 1/2 G. 

Departing from the lF-GLld prescription, we adopt an 
empirical prescription for the energy spectra of the MSPs 
based on the spectra of eight MSPs detected by Fermi, re- 
ported in lAbdo et al.l (|2009al ) . The differential energy spec- 
tra of the Permi-detected MSPs are well-described by a 
power law truncated by an exponential cutoff, 

g«B- r e-^, (6) 

where P is the spectral index and E C ut is the cutoff energy. 
We assume that each spectral parameter, T and E cut , is nor- 
mally distributed in the MSP population, with mean (r) and 
{E cu t) and standard deviation or and o"E cut , respectively. We 
use the spectral parameters of the detected MSPs to iden- 
tify the maximum-likelihood values of these distribution pa- 
rameters, taking into account the measurement uncertainties 
for each pulsar. For this procedure we follo w the methodol- 
ogy described in I Venters Sz Pavlidoul |2007l ). and obtain the 
maximum likelihood parameters [(r),err] = [1.5,0.20] and 
[(_B cut ),cr Bcut ] = [1.9 GeV, 0.54 GeV]. It is notable that the 
distributions are relatively narrow; in particular they imply 
that the vast majority of MSPs have cutoff energies between 
~ 1 and 3 GeV. 

We consider two models for the energy spectra of the 
MSP population. First we examine the simple case in which 
every MSP is assumed to have the same energy spectrum, 
which we denote the Reference Model. In this scenario wc 
assign each MSP the maximum-likelihood average spectral 
parameters, T = 1.5 and E cu t = 1.9 GeV. We also examine 
the impact on our results of allowing the spectral parameters 
to vary within the MSP population according to the distri- 
butions above, and refer to this as the Spectral Variation 
Model. 



3 SIMULATIONS 

Monte Carlo realizations of the Galactic MSP population 
were generated by creating mock MSP catalogs with indi- 
vidual MSP parameters drawn from the distributions given 
in SJ21 To assess the statistical variation between realizations, 
ten Monte Carlo realizations were generated for each case 
considered. The HEALPix package (jGorski et al.ll2005T ) was 
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used to generate maps of the all-sky gamma-ray intensity 
from unresolved MSPs for each mock catalog. Maps were 
constructed at HEALPix order 7 resolution which corre- 
sponds to a pixel size of ~ 0.45° on a side. Each MSP was 
taken to be a point source with no angular extent, and so 
its flux was assigned to a single pixel. 

As we are interested in the emission from unresolved 
MSPs, we excluded from the sky maps the emission from 
MSPs in our mock catalogs which would likely have been 
detected by Fermi. For this purpose we assumed a flux sen- 
sitivity of 10~ 8 p h cm" 2 s" 1 (E > 100 MeV; see, e.g., 
Abd o et al ] |2010al ) and excluded individual MSPs exceed- 
ing this flux threshold. We note that assuming a uniform 
flux sensitivity for MSPs across the sky is a rough approx- 
imation, since point source sensitivity varies with angular 
position due to exposure and foreground contamination, and 
also depends on the individual source spectrum. Although 
we emphasize that this approximation is inadequate for as- 
sessing individual source detectability or completeness, it is 
sufficient for the purpose of removing bright sources which 
are likely to be resolved and would otherwise bias our pre- 
diction for the statistical properties of the diffuse emission. 
Choosing to exclude these sources is conservative, since in- 
cluding bright MSPs would lead to a larger predicted inten- 
sity and a larger contribution to the anisotropy of the IGRB 
from MSPs, and as a result to stronger constraints on MSP 
population models. With this criterion we find that ~ 10 of 
every 10k MSPs in our Reference Model are detectable over 
the entire sky. Note that the parameter JVmsp corresponds 
to the total number of MSPs in the Galaxy and therefore 
includes the detectable sources, although our analysis is per- 
formed on maps of the unresolved sources only. 



4 GAMMA-RAY EMISSION FROM THE 
UNRESOLVED MSP POPULATION 

4.1 Sky distribution of gamma rays from MSPs 

We examined the constraints obtainable on the MSP pop- 
ulation from the intensity and anisotropy properties of the 
Fermi-measured IGRB, so we selected high-latitude sky re- 
gions by excluding Galactic latitudes |6| < 30°. This choice 
matches the latitude mas k applied in the Fermi a ngular 
power spectrum analysis (|Siegal-Gaskins et al. 20101 ). The 
latitude dependence of the results was studied by compar- 
ing the results using a mask excluding |6| < 40°. The choice 
to apply a very generous mask to the Galactic plane also en- 
ables comparison of the high-latitude contribution of MSPs 
to the Fermi-measured IGRB intensity. 

An all-sky map of the gamma-ray intensity for one real- 
ization of the MSP population Reference Model defined in |J2] 
with Amsp = 49k is shown in Figure [T] Emission from in- 
dividual MSPs with fluxes above the detectability threshold 
is not shown. MSPs outside of the latitude mask boundaries 
(marked by lines) are evident, implying a MSP contribution 
to high-latitude diffuse emission. 



4.2 Intensity energy spectra 

The intensity energy spectrum of the gamma-ray emission 
from MSPs outside each latitude mask for the Reference 
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Figure 2. Average intensity energy spectra of the MSP Reference 
Model (solid black line) and Spectral Variation Model (red x's) for 
\b\ > 40° . The intensity spectrum of the Spectral Variation Model 
differs negligibly from that of the Reference Model. The average 
intensity of the MSP Reference Model for |6| > 30° (dashed ma- 
genta line) is also shown. The collective high-latitude intensity of 
the MSPs is more than an order of magnitude smaller than the 
Fermi-measured IGRB intensity (blue crosses) at all energies. 



Model is com pared with the Fer mi-measured IGRB inten- 
sity spectrum l|Abdo et alj|2010b| ) in Fig. [2] The normaliza- 
tion of the MSP intensity outside each mask was obtained by 
averaging over 10 realizations, and the spectral parameters 
of each MSP in the Reference Model were fixed to the maxi- 
mum likelihood values. The average intensity of the emission 
from unmasked MSPs is a factor of ~ 2 larger when exclud- 
ing only \b\ < 30° than when excluding \b\ < 40°, but in 
both cases is more than an order of magnitude smaller than 
the IGRB at all energies. This model is consistent with the 
measured IGRB, but the overall intensity does not provide 
a meaningful constraint on the population. 

In the Reference Model we adopted the simplifying as- 
sumption that all MSPs share the same energy spectrum. To 
test the validity of this assumption, in Fig.[2]we compare the 
collective intensity spectrum of the MSP population, aver- 
aged over |6| > 40°, for the Spectral Variation Model and the 
Reference Model. The differential intensity dN/dE of each 
realization of the Spectral Variation Model was obtained by 
generating a map for each logarithmic energy bin containing 
the integrated intensity of each MSP, given its spectral pa- 
rameters, and then dividing the integrated intensity by the 
energy bin size AE. The points shown represent the average 
intensity outside the mask in each energy bin of 10 Monte 
Carlo realizations of the Spectral Variation Model. To good 
approximation, the collective intensity energy spectrum of 
the Spectral Variation Model matches that of the Reference 
Model, with a small deviation from the Reference Model 
spectrum evident only at the highest energy bin (E ~ 3 
GeV). 

In Fig. [3]we compare the MSP intensity for |6| > 30° to 
the Gal actic diffuse emiss ion for \b\ > 30° (from the model 
used in lCuoco et alJl201Ch . At these latitudes, the intensity 
of the Galactic diffuse emission from cosmic-ray interactions 
with the interstellar gas and photon fields is comparable 
to that of the IGRB, and the MSP emission is subdom- 
inant with respect to both of these signals. However, the 
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Figure 1. MSP gamma-ray intensity integrated from 0.1 to 10 GeV for one realization of the Reference Model. The map is shown 
in Galactic coordinates with the boundaries of the latitude masks excluding |6| < 30° and \b\ < 40° marked. For this figure the map 
resolution was degraded to improve the visibility of MSPs and illustrate their sky distribution; however, all calculations were performed 
on the high-resolution maps as described in the text. 
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Figure 3. Average intensity spectra of the MSP Reference Model 
for |b| > 30° (dashed magenta line) compared with the IGRB 
intensity (blue crosses), t he Galactic diffuse emission for |fe| > 30° 
(dot-dashed red line, from lCuoco et al.l2010l ). and two benchmark 
dark matter models. The two dark matter models correspond to 
an 8 GeV particle pair-annihilating preferentially into t+t - at a 
rate (av) = 1 X 10 _26 cm 3 s _1 (dotted yellow line), and to a 40 
GeV particle annihilating into bb with (av) = 3 X 10 _26 cm 3 s _1 
(solid yellow line). 

Galactic diffuse emission is not expected to contribute sig- 
nificantly to the anisotropy on angular scales of < 1 — 2°, 
corre sponding to multipoles I > 100 (see, e.g.. ICuoco et all 
2010), and therefore MSPs could be a dominant contributor 
to the anisotropy of the high-latitude diffuse emission while 
remaining a subdominant contributor to the intensity. 

Fig. [3] also compares the intensity spectrum of MSPs to 
that of the high-latitude emission predicted for two exam- 
ple dark matter models, chosen because their energy spectra 
bear some resemblance to the collective MSP energy spec- 
trum. We do not resort to any specific particle physics setup 
in the choice of the models. Rather, we specify a dominant 
pair-annihilation final state, the particle mass, and the rate 
of pair-annihilation. One of the dark matter models corre- 



sponds to a dark matter particle with a mass of 8 GeV and 
a cross section {av) = 1 x 10 _26 cm 3 s~ 1 , for which the dom- 
inant annihilation final state is a pair of r leptons. This 
model was chosen to ali g n with that found in the analysis of 
iHooper fc Goodenougbl ||2010D to best fit a gamma-ray ex- 
cess claimed to exist in the innermos t 2 degrees in the direc- 
tion of the Galactic Centre (see also lAbazaiianll2010l for an 
interpretation of that signal as MSP emission). We also com- 
pare a second dark matter model, with a mass of 40 GeV and 
a pair-annihilation cross section (av) = 3 x 10~ 26 cm 3 s _1 , 
for which the dominant annihilation final state is bottom 
quarks. This second model can be regarded as a prototypical 
light bino-like dark matter candidate from the minimal su- 
persymmetric extension of the Standard Model, with a cross 
section that would allow for thermal production of the cor- 
rect universal dark matter density. The intensity of the dark 
matter emission for these two models corresponds to that 
predicted for the high-latitude signal from a nnihilation i n 
Galactic dark matter subhalos in model Al of lAndol {2009), 
assuming the particle properties for each model specified 
above. 

Dark matter annihilation or decay in Galactic substruc- 
ture may generate a significant level of anisotropy in the 
IGRB with an energy dependence similar to that from MSPs 
due to their similar energy spectra. Although the detailed 
shapes of the energy spectra of the dark matter models 
shown in Fig. [3] differ from that of the collective MSP emis- 
sion, the energy range at which both of these possible con- 
tributors become most prominent in the IGRB, as well as 
their cutoff energies, are similar. Since an anisotropy analy- 
sis requires large photon statistics to robustly measure small 
anisotropies, the number of energy bins in which a measure- 
ment can be made with Fermi-LAT is limited, and there- 
fore it may be difficult to localize features in the energy 
dependence of the anisotropy. Consequently, there remains 
the possibility that a MSP-induced anisotropy in the IGRB 
could be confused with a similar signal from dark matter an- 
nihilation. However, we stress that only a signal from very 
light dark matter candidates is likely to exhibit a spectral 
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Figure 4. Angular power spectrum of the Reference Model, with 
Galactic latitudes |6| < 40° masked. Each line corresponds to one 
of ten realizations; the variation between realizations is small. 
The Ci are remarkably constant in multipole, which is consistent 
with the angular power spectrum of an uncorrelated distribution 
of point sources. 

cutoff at sufficiently low energies to effectively mimic MSPs 
in an anisotropy measurement. 

4.3 Angular power spectra 

We consider the angular power spectrum of intensity fluc- 
tuations 5I(ip) = (I(ip) - (/))/(/) where I(ip) is the in- 
tensity in the direction ip, and (7) is the average intensity 
over the unmasked region of the sky. The angular power 
spectrum is calculated by expanding SI in spherical har- 
monics SI = Yli m a i,mYi,m(^) to obtain the coefficients 
| ) . Since a fluctuation map is dimensionless, its 
angular power spectrum characterizes the angular distribu- 
tion of the emission, independent of its overall intensity. 

We calculate the angular power spectrum of the emis- 
sion from MSPs from the simulated sky maps using 
HEALPfx. The angular power spectra are calculated on the 
cut sky, after removing the monopole and dipole compo- 
nents. To approximately correct for the power suppression 
due to masking, the angular power spectra of the cut sky 
are divided by the fraction of the sky outside the mask, 
/ s ky. This approximation is valid at multipoles £ > 100. 

The angular power spectra of 10 Monte Carlo realiza- 
tions of the Reference Model are shown in Fig. [4j calculated 
with a mask excluding jfej < 40°. The scatter between real- 
izations is small, with each realization generating an angu- 
lar power spectrum Ci approximately constant in multipole 
with a value 0.03 < C e < 0.04 for I > 100. The multipole- 
independence of Ci is characteristic of the power spectrum 
of Poisson noise (shot noise) Cp, which arises from an un- 
correlated distribution of sources. Noting that the angular 
power spectrum from MSPs at high latitudes appears to be 
dominated by the Poisson contribution, we hereafter make 
the approximation that the angular power from MSPs is 
constant in multipole, and identify Cp as the average of Ce 
over I = 50 to t = 150. 

The Poisson contribution to the power spectrum scales 
inversely to the number density of sources, i.e., Cp oc 1/Af, 
where A/" is the number of sources per solid angle. Figure [5] 



t 1 1 1 r 




N MSP /1000 

Figure 5. Dependence of Poisson angular power Cp on TVmsp f° r 
the Reference Model, with Galactic latitudes \b\ < 40° masked. 
The error bars represent the standard deviation of the Cp from 
ten realizations, and the solid line illustrates the expected relation 
C P oc l/iV MS p. 

illustrates the dependence of the angular power on the total 
number of MSPs in the model A<msp- Here we adopted the 
Reference Model but varied Nmsp; for each value of A^msp, 
maps for ten realizations of the MSP emission were gen- 
erated. As before, MSPs with individual fluxes above the 
detection threshold were not included in the maps. The av- 
erage Cp of the ten maps for each A^msp is shown in the fig- 
ure. As expected for a Poisson-like source distribution, the 
angular power scales inversely with the number of MSPs. 

It is important to confirm that the expected variation of 
the spectral parameters of MSPs within the Galactic popu- 
lation, in particular the distribution of cutoff energies, does 
not introduce an energy dependence into the angular power 
spectrum. For a source distribution that is independent of 
energy (i.e., a source class in which each member has the 
same observed intensity energy spectrum), the fluctuation 
angular power spectrum is also energy-independent. Energy 
dependence of anisotropy indicates a change in the spatial 
distribution of the contributing sources with energy, and can 
be used to identify the presence of multiple populations or 
populations whose properties var y significantly with energy 
IjSiegal-Gaskins fc Pavlidoull2009r h The angular power spec- 
tra of some astrophysical gamma-ray source populations are 
expected to exhibit a mild, characteristic energy dependence 
due to, e.g., large variations in the spectral properties of indi- 
vidual members of a source class, attenuation by interactions 
with the extraga lactic background light (EBL), and red- 
shifti ng (see, e.g., lAndo fc Komatsu 2006; Z hang fc Beacoml 
l2004t ). In contrast, the angular power spectrum of emission 
from Galactic dark matter annihilation or decay would be 
constant in energy since the energy spectrum is fixed for 
a given dark matter particle, and emission from Galactic 
sources is not subject to redshifting and EBL attenuation. 

The energy dependence of the angular power spectrum 
of the Spectral Variation Model is examined in Fig. [5] The 
angular power spectrum of the MSP emission was calculated 
in each of 10 logarithmically-spaced energy bins, and then 
averaged over ten realizations. The amplitude of the angular 
power spectrum at a given multipole varies by less than a 
factor of 2 over the energy range considered, and the vari- 
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Figure 6. Energy dependence of the angular power spectrum of 
the emission from the MSP Spectral Variation Model. The an- 
gular power spectrum shown for each energy bin has been aver- 
aged over ten Monte Carlo realizations, and was calculated with 
\b\ < 40° masked. 



ation is noticeable only for the highest energy bins (E > 1 
GeV). The slight increase in the angular power at high en- 
ergies is expected for the Spectral Variation Model since, 
due to the variation in cutoff energies within the MSP pop- 
ulation for this scenario, some MSPs no longer contribute to 
the intensity at the highest energy bins, decreasing the num- 
ber density of sources M and thereby increasing the angular 
power Cp. 



5 OBSERVATIONAL CONSTRAINTS 

We now illustrate the potential of anisotropy measure- 
ments to constrain the properties of the MSP population 
by comparing the predictions of our population model to 
the preliminary Fermi anisotropy measurement, and deriv- 
ing constraints on the abundance and emission properties 
of gamma-ray MSPs in the Galaxy. We impose the require- 
ment that the MSPs do not overproduce the IGRB inten- 
sity or anisotropy in the energy range from 1 to 2 GeV, 
and determine the parameter space of MSP models which 
is compatible with this constraint. Since confirmed source 
populations other than MSPs (e.g., blazars) are expected to 
contribute both intensity and anisotropy to the IGRB, al- 
lowing MSPs to contribute all of the measured intensity or 
anisotropy is a conservative choice. 

To assess the contribution of an individual source class 
to the total measured IGRB anisotropy C^tot we construct 
the dimensionful angular power spectrum of the intensity 
by multiplying the fluctuation angular power spectrum of 
a single source class Ce by the mean intensity (7) of that 
source class squared, (I) 2 Ce, where the mean intensity is 
calculated on the unmasked region of the sky. 

Preliminary results from Fermi indicate that the IGRB 
angular power spectrum is approximately constant for £ > 
100, so we identify that measurement as Cp.igrb- The pre- 
liminary Fermi measurement of the IGRB angular power 
spectrum for the energy range of 1-2 GeV, weighted by 
the mean intensity s quared, is (I 2 Cp)igrb — 6.2 x 10~ 18 
(cm -2 s" 1 sr -1 ) 2 sr l|Siegal-Gaskins et al. 201fj| '). The value 



of (I 2 Cp)igrb used here represents the mean of the data 
points from multipoles of £ — 100 to £ = 200. We constrain 
the MSP contribution to the anisotropy at the 2-a level, i.e., 
A 2 ot,MSP C p ,msp ^ (I 2 Cp)igrb + 2o- aniso for 1-2 GeV, where 
Oaniso denotes the mean reported uncertainty on the data 
points. 

Similarly, we derive a constraint from the intensity of 
the IGRB by requiring 7 to t,MSP ^ /igrb + 2cti, integrated 
from 1 to 2 GeV. Fermi's measurement of the IGRB energy 
spectrum is consistent with a power-law with spectral index 
riGR B = 2.41 and I(> lOOMeV) = 1.03 x 10 -5 cm~ 2 s _1 
sr - 1 (|Abdo et al .1 [20 1 bl ^PI. so we adopt this parameteriza- 
tion to determine the integrated intensity of the IGRB from 
1 to 2 GeV, /igrb- The parameter o\ is the reported uncer- 
tainty in the normalization of the power-law fit to the IGRB 
intensity. 

For the MSP intensity and anisotropy, we adopt the 
Reference Model and mask |6| < 30° to match the lati- 
tude mask used in the preliminary Fermi anisotropy mea- 
surement. To explore the dependence of the results on 
the gamma-ray flux distribution of the MSPs, we define 
Fx = /tot, msp ^sk y /(/out Nmsp) to parameterize the typi- 
cal flux contributed by a MSP outside of the mask. The 
parameter /tot.MSP is the mean intensity from all MSPs out- 
side the mask, f2 s k y = 47r/ s ky is the solid angle of the un- 
masked sky, and A^msp is the total number of MSPs in the 
model, as before. The parameter / out is the average frac- 
tion of MSPs outside the mask, which is determined by 
the spatial distribution adopted for the MSPs. We calcu- 
lated /out = 0.03 by averaging the fraction of MSPs with 
\b\ > 30° in 10 Monte Carlo realizations. The intensity 
hot, msp and angular power of intensity fluctuations Cp,msp 
for the MSPs are normalized using the fiducial Reference 
Model values for TVmsp and Fi. The intensity of the MSP 
model then scales as /tot, msp oc FiNmsp and the anisotropy 
as Cp,msp oc I/Wmsp- 

Figure [7] shows the regions of the MSP model pa- 
rameter space excluded by the Fermi measurement of the 
IGRB intensity and preliminary measurement of the IGRB 
anisotropy. The fiducial values of the Reference Model are 
compatible with both the intensity and anisotropy con- 
straints, but an increase in F\ of only a factor of two would 
violate the anisotropy constraint, while the intensity con- 
straint would allow F\ of more than an order of magni- 
tude greater, assuming the fiducial Nmsp- In general, the 
anisotropy constraint is significantly stronger than the in- 
tensity constraint for MSP fluxes and abundances near the 
fiducial values. 

With regard to the intensity constraint, we caution that 
our constraint was derived under the assumption that the 
high-latitude MSP intensity would appear as a contribution 



1 We note that the intensity of the IGRB used in the anisotropy 
measurement is not equivalent to /igrb- The value of /igrb 
was determined by a fitting procedure to remove any spatially- 
dependent components, while the mean intensity of the emission 
for the anisotropy measurement was calculated by simply apply- 
ing a mask excluding \b\ < 30° and masking point sources. The 
procedure we used to calculate the mean intensity of the MSPs, 
i.e., excluding |6| < 30° and emission from individual MSPs above 
the threshold flux, thus closely corresponds to the approach used 
in the Fermi anisotropy analysis. 
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to the me asured IGRB inten sity. The Fermi IGRB mea- 
surement (|Abdo et all |2010bh represents the mean inten- 
sity of the all-sky emission that appears to be isotropic on 
large angular scales. The high-latitude emission from MSPs 
as calculated in the present study is not fully isotropic on 
large angular scales, but rather may exhibit a mild gradient 
away from the Galactic plane. Consequently, although unre- 
solved source populations such as MSPs were not explicitly 
included in the models used in the IGRB intensity spectrum 
analysis, some or all of the unresolved MSP emission could 
have been excluded from the reported measurement of the 
IGRB due to correlations between its large-scale spatial dis- 
tribution and that of the modeled components. 



6 DISCUSSION 

Although tens of thousands of gamma-ray MSPs are ex- 
pected to live in the Galaxy, only a few tens of MSPs have 
now been detected individually at gamma-ray energies, pre- 
senting a significant challenge to constructing an accurate 
population model for this source class. For the purpose of 
studying the potential contribution of Galactic MSPs to the 
anisotropy and intensity of t he IGRB , we adopted the semi- 
empirical models outlined in iF-GLld for the spatial and lu- 
minosity distributions of this population. 

The spatial distribution we used is one that is com- 
monly assumed in semi-analytic population studies of MSPs. 
This spatial model is almost certainly too simplistic, how- 
ever improved versions require population synthesis models 
for Galactic compact objects, including a detailed treatment 
of kinematics in the Galactic gravitational potential and of 
natal supernova kicks, which generally involve severe uncer- 
tainties. Should additional features in the large-scale spatial 
distribution of MSPs be robustly predicted by such studies, 
these could be used as observational signatures tracing a 
MSP component in the diffuse emission. In general, large an- 
gular scale features can provide key information to help iden- 
tify the origin of the observed emission, especially for dis- 
tinguishing a Galactic dark matter signal, which is expected 
to display a spherical symmetry about the Galactic Centre 
on large angular scales, and other Galactic source popu- 
lations, including MSPs, which in stead are typically sym- 
metric about the Galactic plane (|Hooper Sz Serpicq 120071 ; 
llbarra et al.ll2009l ; [Malvshev et al]|20ld ). 

Even without the aid of more detailed models to pin 
down large-scale features in diffuse emission from MSPs, the 
general trend of the dependence of the small-scale anisotropy 
on the spatial distribution of MSPs can be predicted. We 
showed that the Poisson contribution to the angular power 
spectrum Cp is the dominant contribution to the total angu- 
lar power from high-latitude MSPs for multipoles I > 100. 
Since Cp is inversely proportional to the number density 
of sources per solid angle, adopting a model that results 
in a smaller number of MSPs per solid angle in the high- 
latitude sky regions we considered will lead to a correspond- 
ingly larger Cp. 

We demonstrated that, for our Reference Model, the 
fiducial values for the abundance and average high-latitude 
MSP flux imply a non-trivial contribution from MSPs to 
the measured angular power spectrum of the IGRB, despite 
their contribution to the IGRB intensity of only a few per- 



cent. In our treatment we focused on the properties of the 
MSP population which directly impact the IGRB observ- 
ables. We defined the parameter F\ to describe the mean 
flux of a MSP with |6| > 30° in order to encapsulate the 
information contained in the spatial and luminosity distri- 
butions of the MSPs. Casting our results in terms of -Fi 
also provides a means of addressing the uncertainty in the 
gamma-ray efficiency of MSPs, since for MSPs obeying the 
relation L 7 cx E, as the overwhelming majority of MSPs in 
our adopted model do, the gamma-ray flux is linearly pro- 
portional to the assumed efficiency factor /™ ax (see Eq. [3}. 
It is important to note that since the luminosity distribution 
is determined by the distributions of rotation periods and 
magnetic fields within the population, models which predict 
alternative distributions for these parameters would likely 
result in a different predicted Cp. 

To model the energy spectra of the MSP population, we 
used an empirical approach, determining the spectral shape 
and expected distribution of spectral parameters from the 
spectral properties of a small set of Fermi-detected MSPs. 
Based on that sample, we demonstrated that the dimension- 
less angular power spectrum of intensity fluctuations from 
MSPs exhibits only a very mild energy dependence due to 
spectral variation within the population, and that making 
the approximation that all MSPs have identical energy spec- 
tra had a negligible impact on the predicted anisotropy from 
this population. Although the energy-dependence of the di- 
mensionless anisotropy from MSPs alone is minimal, the di- 
mensionful contribution of MSPs to the IGRB anisotropy 
is strongly energy-dependent, since it depends on the rel- 
ative contribution of MSPs to the IGRB intensity at each 
energy. To place constraints, we considered the measured 
IGRB anisotropy at energies from 1 to 2 GeV, which ap- 
proximately maximizes the MSP contribution to the IGRB 
intensity for our assumed MSP energy spectrum. 

A potential challenge for interpreting the results of an 
anisotropy measurement is correctly identifying the source 
populations contributing to the anisotropy. Examining the 
energy dependence of the total anisotropy could help dis- 
tinguish the contributions of different source populations 
by taking advantage of differences in their collective en- 
ergy spectra, but requires that plausible contributing source 
classes have distinct energy spectra. In this sense, MSPs 
could in principle be difficult to distinguish from dark mat- 
ter annihilation or decay in an anisotropy measurement, due 
to a general similarity in the shapes of the MSP and dark 
matter energy spectra for some dark matter models (see 
Fig. [3}- However, the similarity is only likely to be problem- 
atic for very light dark matter candidates with masses < a 
few tens of GeV, and therefore MSPs are unlikely to inter- 
fere with anisotropy-based searches in gamma rays for the 
majority of dark matter candidates. 

In spite of a preliminary Fermi detection of anisotropics 
in the IGRB at energies of a few GeV, the sensitivity of 
anisotropy measurements at higher energies is limited by 
decreasing photon statistics, and consequently constraints 
on anisotropics in the IGRB at the level detected in the 1 to 
2 GeV band are not yet available for energies much above 
this range. We anticipate that future measurements with im- 
proved statistics will be better able to constrain the contri- 
butions of specific source classes to gamma-ray emission us- 
ing energy-dependent anisotropy. Moreover, although we fo- 
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Figure 7. Constraints on MSP population properties from the intensity and anisotropy of the IGRB in the energy range 1-2 GeV. 
MSP models above the lines exceed the measured value of the total IGRB intensity or anisotropy plus 2a. The reference values for the 
parameters Fi and TVmsp are marked. 



cused on the angular power spectrum of the diffuse emission, 
we note that a complementary anisotropy statistic, the 1-pt 
flux PDF, could also be used to help disentangle the con- 
tributions of multiple source populations to diffuse emission 
dLee et alj|2009l ; iDodelson et alj|2009l ; lF-GLld ; iBaxter et al.l 
I2010T ). 



Although we have not considered these signals here, 
anisotropy signatures due to interactions of high-energy 
injected by MSPs may also be imprinted in emission at lower 
energies. This possibility has been pointed out and studied 
in the context of high-energy charged particles produced in 
dark matter annihilation in Galactic subhalos, which could 
lead to signatures in synchrotron emission at radio frequen- 
cies from the intera ctions of in Galactic magnetic fields 
l|Zhang fc Sigll fcoOS). and in low-latitude gamma-ray emis- 
sion from the inverse Compton up-scattering of interstel- 
lar radiation field photons by the injected high-energy e 
l|Zhang et al.ll2010t ). 



As more individual MSPs are detected at gamma-ray 
energies, a more accurate population model will emerge. 
Anisotropy studies can offer complementary information to 
that obtained from studies of detected sources by probing 
the collective properties of the unresolved members of a 
source class. For MSPs, which are a source class with rela- 
tively few individually detected members, anisotropy analy- 
sis could be an important technique to improve our under- 
standing of the characteristics of MSPs at the population 
level, including their abundance and spatial distribution in 
the Galaxy. 
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